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Recent developments in superstring theory and noncommutative geometry are deeply related to the idea of 
Eguchi-Kawai reduction in large N gauge theories which dates back to early 80s. After a general review on this 
subject including revived interests in solving planar QCD, we present some results in the superstring matrix model 
suggesting the dynamical generation of 4d space-time due to the collapse of the eigenvalue distribution. We then 
discuss interesting dynamical properties of held theories in noncommutative geometry, which have been revealed 
by Monte Carlo simulations of twisted reduced models. We conclude with a comment on the recent proposal for 
a lattice construction of supersymmetric gauge theories based on reduced models. 



1. Introduction 

Superstring theory is a unified theory where 
matter and gauge particles are both described by 
various oscillation modes, while all the interac- 
tions including gravity are described simply by 
joining and splitting of strings. Unification of in- 
teractions is not only motivated for aesthetic rea- 
sons but also suggested by precise measurements 
of the coupling constants in the Standard Model. 

In fact superstring theory is so far the only 
quantum theory of gravity that is perturbatively 
well-defined preserving manifest unitarity. This 
is in sharp contrast to the situation with the 
field theoretical approach treating the metric as 
a quantum field, where the theory is perturba- 
tively unrenormalizable. Here one has to study a 
regularized theory nonperturbatively and hope to 
find a nontrivial UV fixed point, where one can 
take the continuum limit. This has been a topic 
studied over a decade. There is still a possibility 
that a sensible continuum limit can be taken, but 
the issue of unitarity remains unclear. 

The problem with superstring theory, on the 
other hand, is that there are too many pertur- 
batively stable vacua with various space-time di- 
mensionality, gauge group, matter contents and 
so on. This means that non-perturbative effects 
are crucial for understanding the 'true vacuum', 
which hopefully describes our world. 



At this point let us recall the history of QCD. 
Properties of its vacuum such as confinement and 
chiral symmetry breaking as well as the dynamics 
of low energy excitations have been understood 
by a nonperturbative formulation of gauge the- 
ory, namely the lattice gauge theory. Likewise 
matrix models, which provide a nonperturbative 
formulation of string theory, are expected to give 
new insights into the nonperturbative dynamics 
of string theory. 

The particular type of matrix models we will be 
discussing appeared in history as large N reduced 
models in the context of solving SU(oo) gauge 
theory £Q . The matrix model proposed as a non- 
perturbative formulation of superstring theory in 
ten dimensions [2] can be regarded as one of such 
models. Among various dynamical issues, we will 
focus on an exciting possibility, which has been 
discussed by many authors [314151617181911011 11121 
113114) . that the SO(10) symmetry of the model is 
spontaneously broken down to SO (4) and the 4d 
space-time appears dynamically. 

String theory has deep connections to noncom- 
mutative geometry, and large N reduced models 
naturally incorporate this feature {15116] . In fact 
a certain type of reduced models j!7j provides 
a lattice regularization of field theories on non- 
commutative geometry [18119120] . We will review 
some Monte Carlo results |21I22I23| . which reveal 
interesting nonperturbative dynamics of such the- 
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ories. 

The reduced model describing nonperturbative 
superstrings has manifest supersymmetry even 
for finite N and this feature has been utilized re- 
cently for a lattice formulation of supersymmetric 
gauge theories |21]. There are also revived inter- 
ests in solving SU(oo) gauge theories with new 
ideas for treating massless fermions and topolo- 
gies of the gauge field |25I26| . In this review we 
also cover these new topics clarifying their mutual 
relationship. 

This article is organized as follows. In Sec- 
tion Owe discuss the connection between matrix 
models and string theory. In Section |2| we in- 
troduce the large N reduced models and discuss 
their equivalence to SU(oo) gauge theories. We 
also review some recent proposals in this direc- 
tion. In Section0]we describe the large N reduced 
model proposed as a nonperturbative formulation 
of superstring theory and discuss in particular the 
dynamical generation of 4d space-time. In Sec- 
tion [5] we discuss the connection between large N 
reduced models and noncommutative geometry. 
We present Monte Carlo results for Yang-Mills 
theory and tfi 4 theory on noncommutative spaces, 
and discuss their intriguing dynamical properties. 
In Section we comment on the relation to the 
recent proposal for a lattice formulation of super- 
symmetric gauge theories. Section |7J is devoted 
to a summary and conclusions. 

2. How matrix models describe strings 

In this Section we review briefly how matrix 
models are related to string theory. For illustra- 
tion let us consider a simple matrix model defined 
by the action 

5 mat = imr0 2 -i^Atr^) 3 , (1) 

where 4> is a N x N Hermitian matrix. This model 
can be solved exactly in the large N limit [27J , but 
in order to see the connection to string theory, we 
make an expansion with respect to the cubic cou- 
pling constant A. Since the dynamical variable 
<t> has two indices, we use the double-line nota- 
tion, which is convenient in discussing the large 
N limit. 




Figure 1. A 3- loop vacuum diagram which ap- 
pears in the perturbative expansion of the (f> 3 ma- 
trix model. 



A typical vacuum diagram which appears as a 
0(A 4 ) contribution to the free energy is shown 
in Fig. ^ Each propagator carries 1/N and each 
vertex carries iVA, according to the Feynman rule, 
which can be read off from the action 0). A 
general vacuum diagram T> with V vertices, P 
propagators and I index loops can therefore be 
evaluated as 

T> ~ fj^j (NX) V N 1 = \ v m . (2) 

An important point here is that the power of N 
is given by the Euler number x = — P + I = 
2(1 — h), which is given by the 'genus' h, namely 
the number of handles of the two-dimensional 
closed orientable surface, on which the diagram 
can be drawn without crossings of lines. For the 
diagram in Fig. we have V = 4, P = 6, I = 4, 
hence h — 0. (Such diagrams are called 'planar 
diagrams' since they can be drawn on a plane 
without crossings of lines as in Fig. Q]) 

As one can see from (01, if one takes the 
large N limit for fixed A, only planar diagrams 
(h = 0) survive. Considering the Feynman dia- 
grams as discretized two-dimensional surface, the 
continuum limit should be taken along with the 
fine-tuning of A ensuring that contributions from 
higher orders in A become dominant. This is 
known to be possible by sending A to the crit- 
ical value A c given by (A c ) 2 = -r^~m GJJ This 
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way of taking the limits (first N — > oo and then 
A y A c ) is relevant to 2d quantum gravity with 
fixed (spherical) topology, if one regards the 2d 
surface as the space-time. If instead one regards 
the 2d surface as the 'worldsheet' describing the 
time-evolution of strings, one may view this the- 
ory as a classical (or tree-level) string theory, 
since there is no branching. 

In this regard let us note that the genus h gives 
the number of loops in the corresponding diagram 
in string theory. In order to obtain all the loop di- 
agrams in string theory, one has to make diagrams 
with any genus survive the large N limit. In the 
present model this is achieved by the so-called 
double scaling limit, which is to take the N — » oo 
limit and the A /" A c limit simultaneously keep- 
ing TV 2 1 A - A c | 5 / 2 fixed I2HESED1- This provides 
a nonperturbative formulation of bosonic (non- 
critical) string theory. In this way one may regard 
matrix models as 'lattice string theory'. 

3. Large N reduced models 

The connection between matrix models and 
string theory that we reviewed in the previous 
Section is quite general. But from now on we will 
focus on a particular class of matrix models which 
is known as large N reduced models. Historically 
they appeared in the context of solving SU(oo) 
gauge theory. 

3.1. Eguchi-Kawai equivalence 

Let us consider SU(iV) lattice gauge theory 1 
on a L D lattice (D-dimensional hypercubic lat- 
tice with the linear extent L) with the periodic 
boundary conditions. Then in the large N (pla- 
nar) limit, Eguchi and Kawai suggested that 
the model with L = oo is equivalent to the model 
with L = 1, namely the one-site model defined by 
the action 

Sek = ~N(3 Y, tr {U^Ulut) . (3) 

1 In this article we will be sloppy about SU(iV) or U(iV) 
to simplify the description since they become equivalent 
in the large N limit. For instance, the in J3J should be 
speciai-unitary, and the symmetry (1) should be {Zj^) D if 
we really start from SU(JV) lattice gauge theory. 



The dynamical variables in this Eguchi-Kawai 
model are given by D unitary matrices of size N . 
(This implies a huge reduction of dynamical de- 
grees of freedom in the large N limit as indicated 
in the title of Ref. p^.) 

In the proof of this statement, there was an 
important assumption that the \](1) D symmetry 

U^e la -U^ (4) 

of the model is not spontaneously broken. 
However it was soon noticed by Bhanot, Heller 
and Neuberger that the U(1) D symmetry is 
actually broken at large (3 (weak coupling) for 
D > 2. This was shown by considering the eigen- 
value distribution of the unitary matrix U^. At 
13 < [3 C the eigenvalues lie uniformly on the unit 
circle, but at (3 > (3 C the eigenvalues are con- 
centrated around some point on the unit circle. 
Since the V(1) D transformation JIJ amounts to 
rotating all the eigenvalues around the unit cir- 
cle, any non-uniform eigenvalue distribution sig- 
nals the SSB. This does not occur at D = 2 in 
accord with the Mermin- Wagner theorem. 

3.2. Remedies to the original proposal 

After the discovery of the U(l) D SSB, reme- 
dies to the original proposal have be suggested. 
The authors of Ref. proposed the quenched 
Eguchi-Kawai model. Their idea was to constrain 
the eigenvalues of to be uniformly distributed 
on the unit circle. This can be achieved by in- 
serting 

J dV^iU^-V^QV}) (5) 

in the partition function, where Q is a diagonal 
unitary matrix Q = diag(l, us, u 2 , • • ■ , w^ 1 ) with 
to = exp(2-7ri / TV) . 

Another proposal was made by Gonzalez- 
Arroyo and Okawa JTj. Their idea was to con- 
sider a L D lattice with twisted boundary condi- 
tions, instead of the periodic ones, and then set 
L = 1. The one-site model thus obtained has the 
action 

5*tek = -N(3 Y V tr {UpUJjfil) , (6) 

where the Z^v factor Z^ v comes from the twist 
in the boundary conditions. This model is called 
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the twisted Eguchi-Kawai model. Modifying the 
boundary conditions does not alter the thermo- 
dynamic limit L — oo, but it does change the 
property of the L = 1 model drastically. The 
configurations which minimize the action are now 
given by which satisfies 

ujfiuF = z^uFupK (7) 

The solution to this equation is unique up to the 
symmetry of the algebra for appropriate choice of 
the twist. Most importantly the minimum-action 
configuration has a uniform eigenvalue dis- 
tribution, which prevents the SSB of U(1) D at 
large (3. The Eguchi-Kawai equivalence of the 
twisted reduced models is reexamined in Ref. [321 • 

Recently Narayanan and Neuberger [JBJ pro- 
posed to consider a partially reduced model with 
L > 1 without quenching or twisting. The key 
observation is that the Eguchi-Kawai equivalence 
holds for arbitrary L as far as V(1) D is not spon- 
taneously broken. The critical /3 = f3 c (L) at 
which the SSB occurs depends on L, and in par- 
ticular (3 C (L) — > oo as L — > oo. Therefore one 
can always choose L such that the (3 one wants 
to study lies below the critical point C (L). 

Various extensions of the Eguchi-Kawai model 
were considered in the 80s. Matter fields in 
the adjoint and fundamental representations have 
been implemented in Refs. 33 34 35 . Extensions 
to non-gauge theories 36 37 38 and to finite tem- 
perature ; 39:;10| are also studied intensively. For a 
comprehensive review on these subjects, we refer 
the reader to Ref. |41| . 

3.3. Revived interests in planar QCD 

In Ref. Kiskis, Narayanan and Neuberger 
proposed to study planar QCD along the idea 
of Eguchi-Kawai reduction. The underlying as- 
sumption is of course that QCD (N c = 3) is ac- 
tually not far from N c = oo as evidenced in the 
glueball mass spectrum. Then the motivation for 
going to the limit N c — oo comes from the fact 
that the valence quark approximation 2 becomes 
exact as far as the number of flavors Nf is kept 
finite. Hence all the pathologies observed with 

2 This is usually called 'the quenched approximation', but 
we prefer not to use this word to avoid confusion with the 
quenching in the reduced model. 



this approximation at A^ c = 3 should go away 
in the N c — ► oo limit, and one can make better 
sense out of Monte Carlo data without having to 
include the dynamical fermions. This is very at- 
tractive given that the full QCD simulation is still 
costly. 

The reason for the exactness of the valence 
quark approximation in the large N c limit can 
be readily seen by comparing the fermion loop 
and the gluon loop. In the double-line notation 
the color index loop running along the gluon loop 
is replaced by the flavor index loop in the case of 
the fermion loop because the fermion is in the fun- 
damental representation with respect to the color 
SU(iV c ). This means that the fermion loop is sup- 
pressed by 0(Nf/N c ) compared with the gluon 
loop. 

If one takes also the Nf — ► oo limit together 
with the iV c — > oo limit with fixed r = Nf/N c , 
which is referred to as the Veneziano limit, the 
effects of sea quarks survive. By varying the pa- 
rameter r, one can smoothly interpolate between 
the valence QCD and the full QCD. 

For the present purpose, one should use the 
quenched Eguchi-Kawai model |31j or the par- 
tially reduced model because the twisted 
Eguchi-Kawai model allows only integer r |35j . 
The authors of Ref. [2E| also suggest to use the 
overlap Dirac operator 021 , which enables the in- 
clusion of exactly massless fermions in the model. 
The problem with the definition for the topologi- 
cal charge mentioned at the end of Section 13.41 is 
also solved by considering the overlap Dirac op- 
erator. Correct chiral anomalies have been repro- 
duced from the reduced models in Refs. |43l44j . 

3.4. Continuum version of reduced models 

Before concluding this Section we comment on 
the 'continuum version' of large N reduced mod- 
els, which is important because manifest super- 
symmetry is implementable. The matrix model 
describing superstrings in ten dimensions falls 
into this category. This version has also been used 
in a lattice construction of supersymmetric gauge 
theories. 

Let us recall that the Eguchi-Kawai model has 
been obtained by considering SU(7V) lattice gauge 
theory and setting the lattice size to L = 1 . Sim- 
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ilarly let us start with continuum SU(7V) gauge 
theory and take the zero- volume limit. The re- 
duced model one obtains in this way is given by 
the action 

S b = ~ix[A li ,A v ] 2 , (8) 

where (fi = 1, • • • , D) are N x N Hermitian 
matrices. 

This model has been studied intensively as the 
bosonic version of the superstring matrix model. 
The finiteness of the partition function (for D > 2 
and sufficiently large N), which is nontrivial be- 
cause the integration region of Hermitian matri- 
ces is non-compact, was shown first numerically 
|45I46| and proved later g7] (see Ref. gH| for a 
recent analytic work on observables). Ref. 05| 
studies the large N behavior of various correla- 
tion functions and in particular the SO(D) sym- 
metry is shown to be unbroken. This means that 
the dynamical generation of space-time expected 
in the D = 10 supersymmetric model does not 
occur in the bosonic version. 

In spite of the importance of such Hermitian 
matrix models in the context of superstrings and 
supersymmetric gauge theories, their equivalence 
to large N gauge theory is dubious. Ref. [50J stud- 
ies the VEV of the Wilson loop numerically in the 
D — 4 bosonic case. The area law holds in a fi- 
nite regime, but this regime neither extends nor 
shrinks in the large N limit. The situation in the 
supersymmetric case is expected to be better |51j . 
but the achieved N (N = 48) is not as large as in 
the bosonic case (N = 768). 

The failure of Eguchi-Kawai equivalence in the 
Hermitian matrix models may be attributed to 
the breaking of the XJ(1) D symmetry, which now 
reads 

A„ i * A M + a M l . (9) 

As in the lattice version, one may consider twist- 
ing or quenching to remedy the situation. It is 
known that twisting is possible only formally at 
N = oo |n2]- Quenching, on the other hand, 
seems to work at least perturbatively, but the 
naive definition of the topological charge Q = 
txF^vF^v vanishes identically \i'6\ . 



4. Lattice superstring 

4.1. The IKKT model 

In 1996 Ishibashi, Kawai, Kitazawa and 
Tsuchiya conjectured that a simple reduced 
model provides a nonperturbative definition of 
type IIB superstring theory in ten dimensions 

0. The model, which is now referred to as the 
IIB matrix model or the IKKT model, can be 
obtained by taking the zero-volume limit of 10- 
dimcnsional SU(iV) super Yang-Mills theory in 
the continuum. The action is given explicitly as 

5'lKKT = Sh + Si , (10) 
Si = — (T.la^WA,,^]) , (11) 

where Sb is the bosonic part JHJ and ty a (a = 

1, • ■ • , 16) are N x N Hermitian fermionic matri- 
ces, which transform as Majorana-Weyl fermions 
under SO(10) transformation. The 16 x 16 matri- 
ces T M are Weyl-projected 7-matrices in the Ma- 
jorana representation. 

There are by now a number of evidences for 
this conjecture. We list some of the most crucial 
ones. Firstly the action l|10l) can be regarded as a 
matrix regularization of the worldsheet action of 
the type IIB superstring theory. Secondly there 
are solitonic objects in string theory, which are 
known as 'D-branes', and the IKKT model con- 
tains these objects as classical solutions. More- 
over the interaction between D-branes calculated 
in the matrix model agrees with the calculations 
in string theory. This is remarkable since the in- 
teraction includes gravity. Thirdly there is an at- 
tempt to derive the string field Hamiltonian from 
Schwinger-Dyson equations in the matrix model 
[53154) . The derivation is successful albeit with 
the aid of a crude power-counting argument. This 
connection, if completed, would provide a direct 
proof of the conjecture. 

There are various dynamical issues that should 
be addressed in this model. Some of them have 
been studied in simplified versions |49l51l55l5fJ] . 
Here we will focus on the possibility that 4d 
space-time appears dynamically in the IKKT 
model. 
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4.2. Emergence of 4d space-time 

Let us first discuss how the space-time is de- 
scribed in this model [3]. For that we diagonalize 
the matrix A,, as 



A» = V^XX 



(12) 



where is a real diagonal matrix X^ = 
diag(xi Al , 352/!, ■ ■ • , £#/»)• By integrating out the 
unitary matrices first, one obtains diagrams 
like Fig. □ but now each index loop is associ- 
ated with Si = {xn, Xi2, • ■ • , Xno), which may be 
viewed as describing the embedding of the world- 
sheet into the 10-dimensional target space. Thus 
we find that the eigenvalues of represent the 
space-time coordinates. This means in particular 
that the space-time is treated dynamically in this 
model. If the eigenvalue distribution of A^ col- 
lapses to a 4d hypersurface in the lOd space-time, 
we are going to obtain 4d space-time dynamically 
in this model. This requires the SO(10) symme- 
try of the model to be spontaneously broken. 

The order parameter for the SSB can be defined 
as follows. Let us define the 'moment of inertia 
tensor' of the space-time as 



T/iv — jytv (A^A V ) 



(13) 



This matrix has 10 real positive eigenvalues, 
which we denote as Ai (i = 1, • ■ • , 10) with the 
fixed ordering 



Ai > A 2 > 



If the ratio 



> Ai 



(14) 



does not approach unity in the 
TV — > oo limit, it signals the SSB. More specifi- 
cally, if (Ai) with i = 1, ■ • ■ , 4 turn out to be much 
larger than (Ai) with i = 5, • • • , 10, it implies the 
dynamical generation of 4d space-time. 

4.2.1. Gaussian expansion method 

In Refs. [10111112] this issue has been addressed 
by the Gaussian expansion method. In general 
the method amounts to considering the action 



Sgem — -^(So 



S - £S ) , 



(15) 



where S is the action of the model one wants to 
study, and So is a Gaussian action. If one sets 
£ = 1 one retrieves the original action S. Then 



one calculates various quantities as an expansion 
with respect to £ up to some finite order and set 
£ = 1. This yields a loop expansion considering 
(So + S) as the 'classical action' and —So as the 
'one-loop counter-term'. The freedom in choosing 
the Gaussian action So is crucial in this method. 
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Figure 2. The free energy of the bosonic model 
obtained by the Gaussian expansion method is 
plotted as a function of v for D = 10. Each curve 
corresponds to the order 1, 2, • • 7. 



Let us first describe how the method works 
in the bosonic version (JSJ of the IKKT model 
|57j . where Monte Carlo results 0H| are available. 
Since we know that the SO(-D) symmetry is not 
spontaneously broken in this model we take 
the Gaussian action to be 

S = ^tv(A fi f (16) 

for simplicity, 3 where the parameter v is left ar- 
bitrary at this point. 

In Fig.JSJwe plot the free energy obtained in the 
Gaussian expansion at various orders as a func- 
tion of the parameter v in (|16f) . Since v is a pa- 
rameter which is introduced by hand, the result 

3 One may also consider SO(-D) breaking Gaussian action 
such as So = v^~^ r (^m) 2 - T ne method works as well, 
and reproduces correctly the absence of SSB in the bosonic 
model. 
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Figure 3. The observable (■^tr(A /x ) 2 ) in the 
bosonic model obtained by the Gaussian expan- 
sion method at orders 1, • • • , 7 for various D. The 
horizontal lines represent either the Monte Carlo 
results (D < oo) or the exact result (D = oo). 



should not depend much on it if the expansion be- 
comes valid for some value of v. Remarkably we 
do see a clear plateau in Fig. EI and the height of 
the plateau stabilizes with the increasing order. 
Similar behaviors were obtained for the observ- 
able (-^tr (A M ) 2 }, and the height of the plateau 
estimated at each order (with the 'error bar' rep- 
resenting the uncertainty) is shown in Fig. EI The 
results converge to the horizontal lines, which rep- 
resent either the Monte Carlo results (D < oo) or 
the exact result (D — oo). The convergence be- 
comes faster for larger D. 

Let us turn to the IKKT model (jTU|) . and de- 
scribe how the first evidence for the 4d space- 
time was obtained in Ref. ^Oj- Since the SSB of 
the SO(10) symmetry is the main issue here, the 
Gaussian action 

s o = £ — tr [A^f + NA aP tx (* Q ^) , (17) 

p V » a/3 

which breaks the SO(10) invariance, was consid- 
ered, and the Gaussian expansion has been per- 
formed up to the 3rd order. Since it is difficult to 
search for plateaus in the whole space of and 
A a p, it was assumed that SO(ri) times some dis- 



Figure 4. The free energy obtained by the Gaus- 
sian expansion method for the IKKT model at 
orders 1,3 and 5. The two symbols correspond 
to the SO(4) Ansatz and the SO(7) Ansatz re- 
spectively. At the 3rd order and higher the free 
energy for the SO (4) Ansatz becomes lower than 
the other. 



crete subgroup of SO(10 — d), where d = 2, 4, 6, 7, 
is preserved. This assumption reduced the num- 
ber of parameters to three, and plateaus were 
searched for by finding extrema of the free en- 
ergy with respect to the three parameters for each 
d. At order 3, it was found that the solution for 
d = 4 gives the smallest free energy. The ex- 
tent of the space-time in the d dimensions (R) 
and that in the remaining (10 — <i) dimensions (r) 
have been estimated at the points in the param- 
eter space which extremize the free energy. The 
ratio R/r turned out to be larger than one, and 
for d = 4 it increased drastically as one goes from 
order 1 to order 3. 

A lot of effort has been made to increase the 
order of the expansion. It was soon noticed that 
Schwinger-Dyson equations can be used to reduce 
the number of Feynman diagrams to be evaluated 
|llj . Then a computer code has been written in 
order to automatize the task of listing up and 
evaluating all the Feynman diagrams ^2j- With 
these technical developments, the order of the ex- 
pansion has now been increased up to the 7th or- 
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has been simulated in Ref. 0]. We shall denote 
the VEV with respect to this partition function 
by ( • }o- The results for the order parameters 
(Ai)o are shown in Fig-El If one makes a linear ex- 
trapolation to N — oo, one finds that all the (Aj)o 
converge to the same value, meaning that there 
is no SSB in the phase-quenched model. This im- 
plies that the phase of the fermion determinant 
plays a crucial role in the SSB if it happens at all 
(Refs. |5I8| support this scenario). 



0.035 




Figure 5. The ratio R/r of the extents obtained 
by the Gaussian expansion method for the IKKT 
model at orders 1,3 and 5. The two symbols corre- 
spond to the SO(4) Ansatz and the SO(7) Ansatz 
respectively. The result for the SO (4) Ansatz 
grows with the order. 



der, and the results strengthened the conclusion 
of Ref. ^01- I n Figs. 0] and [3] we show the free 
energy and the ratio R/r of the extents obtained 
up to the 5th order 1 1 . 

4.2.2. Monte Carlo simulation 

Although the results of the Gaussian expan- 
sion method are encouraging and deserve further 
investigations, it is also important to confirm the 
results by Monte Carlo simulation from first prin- 
ciples. This approach will also allow us to un- 
derstand the mechanism for the collapsing of the 
eigenvalue distribution of . 

An important point to note here is that the 
fermion determinant det M — | det M \ e lT is actu- 
ally complex because the fermionic matrices \& a 
transform as lOd Major ana- Weyl fermion, which 
is essentially chiral. (In Euclidean space chiral 
determinants are generally complex.) This poses 
the notorious technical problem known as the 
'complex-action problem'. As a first step, the 
phase-quenched model 
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Figure 6. The 10 eigenvalues of the moment of 
inertia tensor in the phase- quenched model are 
plotted against 1/N. 



In order to study the effect of the phase, a new 
method has been considered HJ. Let us first nor- 
malize the eigenvalues of the moment of inertia 
tensor as 



Ai = 



(19) 



Z 



dAe 



det Ml 



(18) 



(Ai)o ' 

Then we define the distribution 

Pi(x) = (6(x- Ai)) , (20) 

which has the factorization property 

Pi(x) = ±p?\x)w i (x) , (21) 

where pf\x) = {5(x — A,))o denotes the distri- 
bution defined in the phase-quenched model (|18l) 



and C = (e l 



is a normalization constant. 
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The weight factor Wi(x), which represents the 
effect of the phase, can be obtained by performing 
the constrained simulation 



Zi(x)= I dAe~ Sb \detM\6(x- X l ) 



(22) 



and calculating the expectation value of e lF . This 
calculation is the most time-consuming part be- 
cause of the oscillating phase, but one can still 
do it for moderate N. Various virtues of the 
method, as compared with the standard reweight- 
ing method, are discussed in Ref. |55| . 



o 
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Figure 7 

plotted for various N. The solid curve represents 
a fit to some analytic function. The dotted curve 
represents a similar fit to the data for $5 (a;). 



The distribution pf^ (x) for the phase-quenched 
model is peaked at x — 1 due to the chosen 
normalization Ijl9|l . The weight factor Wi(x), on 
the other hand, turned out to have a minimum 
around x ~ 1 (see Fig. [7J), and as a result the 
distribution pi (x) for the full model has a double- 
peak structure. Which of the two peaks becomes 
dominant at large N may depend on i. If it turns 
out that the peak at x > 1 dominates for i < 4, 
but the peak at x < 1 dominates for i > 5, we are 
able to obtain 4d space-time. 
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(23) 



A crucial point in this approach is that 
1 

scales with N as one can see from Fig. [7| This 
scaling behavior is also understandable by gen- 
eral arguments |^]- Using the scaling function 
<&i(x) extracted in this way, one can estimate 
the height of the two peaks in the distribution 
Pi(x) at much larger N. Preliminary results for 
N = 64, 128 are encouraging, but it remains to 
be seen whether one can definitely conclude that 
4d space-time appears in the IKKT model. 

5. Lattice noncommutative geometry 

In this Section we discuss noncommutative ge- 
ometry, which has deep connections to string the- 
ory. In particular it was shown by Seiberg and 
Witten jSH] that field theories in noncommuta- 
tive geometry appear as a certain low-energy limit 
of string theory, and this triggered a tremendous 
boom in this subject. Among other things it was 
realized that such theories have various dynami- 
cal properties that ordinary field theories do not 
have. This is due to the UV/IR mixing effect, 
which was discovered in the one-loop calculation 
|60| . Calculations beyond one-loop become com- 
plicated because of this effect, and perturbative 
renormalizability is not proved even for simple 
scalar field theories [HI]- This raises some sus- 
picion that these theories are actually not well- 
defined nonperturbatively or the UV/IR mixing 
effect is merely a one-loop artifact. 

Fortunately now we know how to regularize 
these theories on the lattice. In fact it was found 
that twisted reduced models at finite N can be 
interpreted as lattice-regularized field theories in 
noncommutative geometry |19j . This is a refine- 
ment of the earlier work |lfi| . which pointed out 
the connection between twisted reduced models 
and noncommutative field theories. 

Noncommutative geometry is characterized by 
the commutation relation among the space-time 
coordinates 



where Hermitian operator. 



(24) 

When we re- 



place ordinary coordinates x^ by the noncommut- 



10 



ing ones x M , the field cj>(x) should be replaced by 
$ = 4*{x), which is also an operator. 

In fact this setting can be put on a periodic 
L D lattice in such a way that there is a one-to- 
one correspondence between a field <fi(x) on the 
lattice and a N x N matrix <!>. This means in 
particular that 



L D =N 2 



(25) 



from the matching of the dynamical degrees of 
freedom. Using this matrix-field correspondence, 
one can derive the noncommutative lattice field 
theories from twisted reduced models [TJJ. As- 
suming the canonical form for the noncommuta- 
tive tensor 9^ v in l(2*I)l. its scale 4 is found to be 



-La 2 



(26) 



where a is the lattice spacing. 

If one takes the planar limit sending N to in- 
finity with fixed a, one obtains 9 = oo. On the 
other hand, the twisted reduced model becomes 
equivalent to large N field theories in this limit 
as we discussed in Section 13.21 This is a non- 
perturbative account of the well-known fact that 
noncommutative theories at 9 = oo are equivalent 
to large N theories, which can be easily shown di- 
agrammatically (see footnote [SJ . 

In order to obtain noncommutative theories at 
finite 9, one has to take the large N limit and 
the continuum limit simultaneously. More specif- 
ically one has to take the limits N — > oo and 
a — > in such a way that N 2 / D a 2 is fixed. This 
corresponds to the double scaling limit we dis- 
cussed in Section [21 In what follows we discuss 
the dynamics of specific theories, which have been 
studied by Monte Carlo simulation. 

5.1. NC Yang-Mills theory on the lattice 

As the simplest possible gauge theory in non- 
commutative geometry, 2d pure Yang-Mills the- 
ory has been studied by Monte Carlo simulation 
[2*T] using the twisted Eguchi-Kawai model J§J|. 
The planar limit of the theory is solved on the 
lattice by Gross and Witten |£2|, so we can refer 
to the result in our analysis. Let us consider the 

4 In the D = 2 case, for instance, is defined by M „ = 



VEV of the Wilson loop (W(IxI)), which is com- 
plex in general unlike in ordinary gauge theory 
because the noncommutativity 9 breaks parity. 
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Figure 8. The polar coordinates of the complex 
Wilson loop W(I) plotted against the physical 
area A = a 2 1 2 . At small areas it follows the 
Gross- Witten area law (solid line). At larger ar- 
eas the absolute value does not decay any more, 
but the phase increases linearly. 

In Fig. ISI we plot the absolute value and the 
phase of the Wilson loop obtained for fixed 9 = — . 
One can see that the data in the small area regime 
(al) 2 <§; 9 agree with the Gross- Witten planar re- 
sult. This is understandable because in the UV 
regime one cannot see the finiteness of 9, so the 
result becomes indistinguishable from that with 
9 = oo, namely the planar result. The coupling 
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constant (3 has been fine-tuned as a function of 
the lattice spacing a in such a way that the re- 
sults in the planar regime scale. (In the present 
case we can use the Gross- Witten result to infer 
the tuning of the coupling constant.) The scaling 
in this regime is simply a consequence of the fact 
that the commutative large N gauge theory has a 
sensible continuum limit. What is highly nontriv- 
ial is that the scaling extends to larger area, where 
nonplanar diagrams start to contribute. This rep- 
resents the continuum limit of the noncommuta- 
tive theory with finite 9. 

In particular we find that the phase of the Wil- 
son loop grows linearly with the area as 



(phase) = 9 1 x (area) , 



(27) 



which is reminiscent of the Aharonov-Bohm ef- 
fect if we identify 9~ 1 with a static magnetic field 
traversing the 2d plane. Such an identification 
occurs commonly in noncommutative geometry, 
but the dynamical AB effect observed here awaits 
more profound understanding. Let us also com- 
ment that the IR behavior represented by 1)27(1 is 
clearly different from ordinary gauge theory. This 
confirms that the introduction of the noncommu- 
tativity 9 can change the IR physics, which is not 
expected at the classical level. One should un- 
derstand this fact as a consequence of the UV/IR 
mixing effect caused by the nonplanar diagrams. 

5.2. NC cj) A theory on the lattice 

Let us discuss the effects of nonplanar diagrams 
more closely in the case of scalar field theories. 
For the one-loop correction to the inverse propa- 
gator, one obtains 



F(p) = d?q 



q 2 + to 2 



(28) 



where the unusual phase factor e t9 ^^p^ [ s ^ ne 
only effect of the noncommutative geometry. 5 

Because of the 'noncommutative phase', the in- 
tegration over the loop momentum q converges, 



5 In the case of planar diagrams, this 'noncommutative 
phase' cancels. In the 9 — ► oo limit, all the nonplanar 
diagrams vanish due to the oscillating phase, and only 
the planar diagrams survive. This is the diagrammatic 
account for the equivalence of the 9 = oo theory to the 
large N field theory. 
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Figure 9. The phase diagram of the (2 + l)d 
noncommutative <j> theory on the lattice. 



and the function F(p) becomes finite for p =/= 0. 
At p — the phase vanishes, and the momentum 
integration has the usual UV divergence. This is 
reflected in the singularity of F(p) at p ~ as 



F(p) 



1 



\9p\ 



D-2 



(29) 



Thus the UV divergence of the commutative the- 
ory is transformed into the IR singularity with 
respect to the external momentum. This is the 
famous UV/IR mixing effect 60 . 

As a physical consequence of this effect, Gubser 
and Sondhi j^SI conjectured the existence of the 
striped phase, where a nonzero momentum mode 
condenses, thus giving a position dependent VEV 
to the scalar field. This in particular means that 
the translational invariance is spontaneously bro- 
ken. The analysis was based on self-consistent 
Hartree-Fock approximation at the one-loop, and 
it would be interesting to examine this conjecture 
by Monte Carlo simulation. 

Let us consider (2+l)d NC </> 4 theory on the 
lattice |22j . where two spatial directions are non- 
commutative as in Section 15.11 but now we also 
have a commuting (Euclidean) time direction in 
addition. The matrix model describing such a 




Figure 10. The dispersion relation in the disor- 
dered phase near the critical point. 



theory is given by 



Ntr 



~X;{iV*(t)rt-*(t)} ! 



/x=l 



+ -($(£ + l)-<I>(i)) 2 



(30) 



where the shift operators in the noncommut- 
ing directions are defined by 



r r 

z, u . 



■ Af + l 

e 7 ™— = Z* 



(31) 
(32) 



Fig. El shows the phase diagram obtained by 
Monte Carlo simulation. The disordered phase 
appears at large m 2 as in the well-known commu- 
tative case. Let us then decrease m 2 for a fixed 
coupling constant A. When A is small, the system 
undergoes the usual Ising-type phase transition 
into the uniformly ordered phase. However, when 
A is sufficiently large, one obtains the striped 
phase as conjectured in Ref. Whether the 

width of the stripes becomes finite in the contin- 
uum limit is an interesting open question, which 
is currently being investigated. 

This phenomenon can be understood more 
clearly by looking at the dispersion relation. Fig. 
1101 shows the dispersion relation obtained in the 



Figure 11. A snapshot of the configuration ob- 
tained just below the critical point with the same 
A as in Fig. The average has been taken over 
the X2 direction, in which the value of the scalar 
field is uniform. 



disordered phase near the critical point. One can 
see an anomalous behavior at p ~ 0, and as a re- 
sult the minimum of the energy occurs at nonzero 
momentum. The data can be fitted to the form 



E 2 =p 2 



M 



A 

W\ 



(33) 



suggested by the one-loop calculation Fig. 
II H is a snapshot of the configuration obtained just 
below the critical point. The data can be nicely 
fitted to the sine function, in accord with the con- 
densation of the nonzero momentum mode. If m 2 
is even lower, the shape becomes more step-like. 

6. Lattice supersymmetry via orbifolding 

In this Section we comment on the recent pro- 
posal for a lattice construction of supersymmetric 
gauge theories based on reduced models . Let 
us recall that the continuum version of reduced 
models can incorporate manifest SUSY even at 
finite N. Here the anti-commutator of Q and Q 
yields a generator of the U(l) D transformation 
(JSJ. Thus the translation is realized in the inter- 
nal space, which is consistent with the identifica- 
tion of the eigenvalues of as the space-time 
coordinates described in Section FQ1 

The key to get a lattice theory from the con- 
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tinuum reduced model is the 'orbifolding', which 
is to impose a constraint on a matrix $ such as 

n^n] = e ir > $ for j = 1, ■ • • , d . (34) 

Solving the constraint explicitly, one obtains a 
d-dimensional lattice field theory. In particular 
the 'charge vector' rj in (|34|l determines that the 
lattice field corresponding to the matrix $ lives 
on the links connecting x and x + r . 

Part of the SUSY survives the orbifolding, 
which enables the restoration of the full SUSY in 
the continuum limit without fine-tuning. An ex- 
tension to the noncommutative geometry is pos- 
sible by making Clj in 134(1 noncommutative |64| . 

7. Summary 

We have seen that the idea of Eguchi-Kawai 
reduction in large N gauge theories has developed 
in many different directions. Let us summarize 
them classifying large N reduced models into the 
lattice version and the continuum version. 

In the lattice version there are a few ways to 
achieve the equivalence to large N gauge theo- 
ries, and there are revived interests in this direc- 
tion using the quenching or the partial reduction. 
On the other hand, the twisted reduced mod- 
els have been given a new interpretation at finite 
TV as a lattice regularization of noncommutative 
field theories. This motivated a different type of 
large N limit other than the planar limit, where 
non-planar diagrams also survive. Monte Carlo 
studies revealed intriguing dynamical properties 
of these theories caused by the UV/IR mixing. 

The continuum version, on the other hand, has 
been obscure in the context of equivalence to large 
N gauge theories, but it suddenly became im- 
portant as a nonperturbative formulation of su- 
pcrstring theory. We have discussed that the 
collapsing of the eigenvalue distribution, which 
was problematic in the context of Eguchi-Kawai 
equivalence, may provide the key to understand 
the dynamical generation of 4d space-time in lOd 
superstring theory. The ultimate goal in this di- 
rection is of course to derive the Standard Model 
from the IKKT model or from whatever matrix 
models describing nonperturbative superstrings. 
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